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Abstract 

We present improved worldline numerical algorithms for high-precision calculations of 
Casimir interaction energies induced by scalar-field fluctuations with Dirichlet bound- 
ary conditions for various Casimir geometries. Significant reduction of numerical cost 
is gained by exploiting the symmetries of the worldline ensemble in combination with 
those of the configurations. This facilitates high-precision calculations on standard 
PCs or small clusters. We illustrate our strategies using the experimentally most 
relevant sphere-plate and cylinder-plate configuration. We compute Casimir curva- 
ture effects for a wide parameter range, revealing the tight validity bounds of the 
commonly used proximity force approximation (PFA). We conclude that data anal- 
ysis of future experiments aiming at a precision of 0.1% must no longer be based on 
the PFA. Revisiting the parallel-plate configuration, we find a mapping between the 
D-dimensional Casimir energy and properties of a random-chain polymer ensemble. 

1 Introduction 

Recent years have witnessed remarkable qualitative and quantitative progress in the un- 
derstanding of the Casimir effect P both experimentally and theoretically. Measurements 
of the Casimir force have reached a precision level of 1% El HI El El U\- Further im- 
provements are currently aimed at with intense efforts, owing to the increasing relevance 
of these quantum forces for nano- and micro-scale mechanical systems. Also from the per- 
spective of fundamental physics, Casimir precision measurements play a major role in the 
search for new sub-millimeter forces, resulting in important constraints for new physics 

[ainiiiniiniini. 

On this level of precision, corrections owing to surface roughness, finite conductivity, 
thermal fluctuations and geometry dependencies have to be accounted for (THl EU EE 
El El- These corrections may be classified in terms of material corrections on the one 
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hand; they are induced, for instance, by surface roughness and finite conductivity which 
may be viewed as a deviation from the ideal Casimir configuration. On the other hand, 
corrections due to geometry dependence are of direct quantum origin and thus universal, 
i.e., independent of the microscopic details of the interactions between the fluctuating field 
and the constituents of the surfaces. Since material corrections are difficult to control with 
high precision, force measurements at larger surface separations up to the micron range 
are intended. 1 Though this implies stronger geometry dependence, this universal effect is, 
in principle, under clean theoretical control, since it follows directly from quantum field 
theory [H]. 

Straightforward computations of geometry dependencies have long been conceptually 
complicated, since the relevant information is subtly encoded in the fluctuation spectrum. 
Generically, analytic solutions are restricted to highly symmetric geometries. This problem 
is particularly prominent, since current and future precision measurements predominantly 
rely on configurations involving curved surfaces, such as a sphere above a plate. Curved 
configurations help to circumvent the difficulty of maintaining parallelism as it occurs in 
the parallel-plate configuration; the latter has been mastered so far only in one experiment 
[22] with a precision level of ~ 15%. As a general recipe for curved configurations, the 
proximity force approximation (PFA) [22] has been the standard, though uncontrolled, tool 
for estimating curvature effects for non-planar geometries in all experiments so far. 

In recent years, various new techniques have been developed for computing Casimir 
effects in more involved geometries [2111231211123121112312311]], each with its own 
merits and limitations. This includes improved approximation techniques which can deal 
with curved geometries more reliably, such as the semiclassical approximation [21] and the 
optical approximation [21], as well as exact field theoretic methods based on functional- 
integral techniques [23 123 123 or scattering theory [23 12H] • 

In this work, we use and further develop worldline numerics [23 , which facilitates 
Casimir computations from field-theoretic first principles. Worldline numerics builds on a 
combination of the string-inspired approach to quantum field theory [22] with Monte Carlo 
methods. As a main advantage, the worldline algorithm can be formulated for arbitrary 
geometries, resulting in a numerical estimate of the exact answer [23- The inherent use 
of Feynman path-integral techniques circumvents the problem of determining the Casimir 
fluctuation spectrum [32], which is often encountered in other approaches. The resulting 
algorithms are trivially scalable and computational efforts increases only linearly with the 
parameters of the numerics. 

Here, we present worldline algorithms to examine the Casimir effect in a sphere-plate 
and cylinder-plate geometry for a fluctuating scalar field, obeying Dirichlet boundary con- 
ditions ("Dirichlet scalar"). We compute the Casimir interaction energies that give rise to 
forces between the rigid surfaces. This allows for a quantitative determination of validity 
bounds of approximation methods such as the PFA, some results of which have already 

1 Measurements at larger surface separations are also aimed at in order to resolve a recent controversy 
about thermal corrections, see |18U19j and references therein. Even though thermal contributions are also 
universal in the ideal Casimir limit, they can mix nontrivially with material corrections in a way that may 
affect any real experiment |2()j . 
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been presented in a recent Letter [31]. We detail significant improvements of the numerical 
algorithms which facilitate high-precision calculations. Apart from numerical discretization 
errors which are kept at or below the 0.1% level, no quantum-field-theoretic approximation 
is needed. Our results further strengthen the agreement with recently obtained analytic 
solutions for medium or larger curvature [2H1 EHj and for small curvature [HD] • 

We emphasize that the Casimir energies for the Dirichlet scalar should generally not 
be taken as an estimate for those for the electromagnetic (EM) field, leaving especially 
the experimentally most used sphere-plate case as a pressing open problem. Neverthe- 
less, a comparison with other techniques can meaningfully be performed, and the validity 
constraints that we derive, e.g., for the PFA hold independently of the type of boundary 
condition, since the PFA approach makes no reference to the nature of the fluctuating 
field. If an experiment is performed outside the PFA validity ranges determined below, 
any comparison of the data with theory using the PFA has no firm basis. 

We also revisit Casimir's classic parallel-plate configuration, first because further algo- 
rithmic strategies can easily be illustrated here; and second, we thereby obtain a mapping 
between the D dimensional Casimir effect and characteristic properties of a random-chain 
polymer ensemble (say, in 3 space dimensions), due to the use of path integrals in the 
worldline method. 

Our work is organized as follows: in Sect. El we briefly review elements of the worldline 
formulation for Casimir configurations as well as the basic ideas of worldline numerics. In 
Section El the construction of our new worldline algorithms is detailed for various Casimir 
configurations. We present our conclusions in Sect. 0] We close this introduction with 
a brief review of Casimir curvature effects and the PFA; the latter is not only a simple 
(though potentially misleading) approximation, but also provides for a useful normalization 
for our numerical results. 

1.1 Casimir curvature effects and proximity force approximation (PFA) 

An intriguing property of the Casimir effect has always been its geometry dependence. As 
long as the typical curvature radii Ri of the surfaces are large compared to the surface 
separation a, the PFA is assumed to provide for a good approximation. In this approach, 
the curved surfaces are viewed as a superposition of infinitesimal parallel plates [23*1 116j. 
The Casimir interaction energy is obtained by an integration of the parallel-plate energy 
applied to the infinitesimal elements. Part of the curvature effect is introduced by the 
choice of a suitable integration measure which is generally ambiguous, as discussed, e.g., 
in j27j. For the case of a sphere with radius R at a (minimal) distance a from a plate, the 
PFA result at next-to-leading order reads 




(1) 
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where the upper (lower) coefficient in braces holds for the so-called plate-based (sphere- 
based) PEA. They represent two limiting cases of the PFA and have often been assumed to 
span the error bars for the true result. Furthermore, cpp = 2 for an EM field or a complex 
scalar, and cpp = 1 for a real scalar field. 

The first field-theoretic confirmation of the zeroth-order result £pp A (a, R) has been 
obtained within the semi-classical approximation in [21]. We will use this zeroth-order 
interaction energy (and its analogue for the cylinder-plate configuration) as a normalizer 
for our numerical estimates. As an advantage, any deviation from this result can be 
interpreted as a true quantum- induced Casimir curvature effect. Future experiments are 
indeed expected to become sensitive to the first-order curvature correction, which therefore 
is of particular interest to us. 

Conceptually, the PFA is in contradiction with Heisenberg's uncertainty principle, since 
the quantum fluctuations are assumed to probe the surfaces only locally at each infinites- 
imal element. However, fluctuations are not localizable, but at least probe the surface in 
a whole neighborhood. In this manner, the curvature information enters the fluctuation 
spectrum. This quantum mechanism is immediately visible in the worldline formulation of 
the Casimir problem. Therein, the sum over fluctuations is mapped onto a Feynman path 
integral, see below. Each path (worldline) can be viewed as a random spacetime trajectory 
of a quantum fluctuation. Owing to a generic spatial extent of the worldlines, the path 
integral directly samples the curvature properties of the surfaces [2E]. 

2 Casimir effect on the worldline 

2.1 Worldline formulation for a Dirichlet scalar 

Let us briefly recall from j2B] how the Casimir effect for a real scalar field <fi satisfying 
Dirichlet boundary conditions can be described in the worldline formalism. For this, the 
field is coupled to a background potential V (x) > which models the Casimir configuration: 
qualitatively, the amplitude of the field fluctuations are suppressed at those regions of 
spacetime where the potential V(x) is large. The field-theoretic Euclidean Lagrangian is 



Here, we have included a mass term for the scalar; the massless limit, which is more 
analogous to the photon field, can always safely be taken. From the standard generating 
functional for the quantum correlation functions of V(x), we obtain the (unrenormalized) 
effective action, which reads in worldline representation in D = d+1 dimensional Euclidean 
spacetime: 
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The expectation value (. . . ) x has to be taken with respect to an ensemble of closed world- 
lines with GauBian velocity distribution, 



(...)*:= / V.r t J"''- ( 5 ) 



x(T)=x(0), 
CM 



with implicit normalization (l) x = 1. For convenience, the common center of mass xcm of 
the worldlines has been separated off. 

In this work, we concentrate on Casimir forces between disconnected rigid bodies which 
we represent by a time-independent potential V(x) = Vi(x) + V 2 (x) + ■ ■ ■ ; the potentials 
Vi(x) for the single bodies have pairwise disjoint supports, i.e., ^(x)V^(x) = for all x and 
% j. From the effective action, we obtain the Casimir energy by scaling out the trivial 
time integration, 

£ = TT-- ^ 

J «^0,CM 

For the Casimir force, only the portion of the Casimir energy which depends on the relative 
position of the objects is relevant. This portion can conveniently be extracted from the 
total Casimir energy by subtracting the (self-)energies of the single objects. This leads us 
to the Casimir interaction energy, 

-^Casimir := Sv 1 +V 2 +- ~ &Vi ~ ^V 2 ~ ' " " > (?) 

which serves as the potential energy for the Casimir force; i.e., Casimir forces (or torques, 
etc.) are obtained by the (negative) derivative of -Ecasmur with respect to a distance (or 
angle) parameter. By this procedure, also any UV divergencies of Eq. (0) are automatically 
removed. Moreover, the interaction energy can thus be well defined even if the Casimir 
(self-)energy of a single surface is ill-defined in the ideal boundary-condition limit ("perfect 
conductivity", infinitely thin surfaces, etc.) j35J EHl EZj • 

For the ideal limit of infinitely thin surfaces, the potential V(x) becomes a 5 function 
in space, 

V(x) = xj da 5 (d) (x-x CT ). (8) 

The geometry of the Casimir configuration is defined by E, denoting a d — 1 dimensional 
surface. The surface measure da is assumed to be re-parameterization invariant, and x^ 
denotes a vector pointing onto the surface. For a typical configuration, S consists of two 
disconnected objects (e.g., two disconnected plates), £ = Si U S 2 , with Si PI S 2 = 0. The 
positive coupling A has mass dimension 1. In the ideal limit A — > oo, the potential imposes 
Dirichlet boundary conditions on the quantum field. 

For the potential Eq. (jSJ), the r integral in the expectation value in Eq. (J3J) reads 

Wr)] := [ T dr V(x(r)) = A £ — ±— , (9) 
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where the sum goes over all intersection points of the worldline x(T a ) and the surface £. 
In the denominator, xi(r ff ) denotes the component of the r derivative perpendicular to 
the surface. 

Computing the Casimir interaction energy Eq. (JZj) for two surfaces Si and £2, the 
argument of the expectation value in (J3J becomes 

( e -/E lU E 2 [x(r)] _ y _ ( e -J El [x(r)] _ ^ _ ( e -/E 2 [x(r)] _ ^ e [ Q) ^ ( 1Q ) 

Most importantly, Eq. is nonzero only if the loop x(r) intersects both surfaces. In 
the Dirichlet limit A — * 00, this expression then equals one. Thus, for a massless scalar 
field with Dirichlet boundaries in D = 3 + 1, the worldline representation of the Casimir 
interaction energy boils down to IHH] 

11 f°° dT f 

^Casimir = "3 (4^2 ] J*] ^ * CM ^W^x' ( U ) 

Here, the worldline functional Gs[x(r)] = 1 if the path x(r) intersects the surface S = 
Si U S 2 in both parts Si and S 2 , and 6s[x(r)] = otherwise, analogous to the standard 
step function. 

This compact formula has an intuitive interpretation: the worldlines can be viewed as 
the spacetime trajectories of the quantum fluctuations of the <fi field. Any worldline that 
intersects the surfaces does not satisfy Dirichlet boundary conditions. All worldlines that 
intersect both surfaces thus should be removed from the ensemble of allowed fluctuations, 
thereby contributing to the negative Casimir interaction energy. The auxiliary integration 
parameter T, the so-called propertime, effectively governs the extent of a worldline in 
spacetime. Large T correspond to IR fluctuations with large worldlines, small T to UV 
fluctuations. Those T values for which the spatial extent of the worldlines is just big enough 
to intersect with both surfaces generically dominate the Casimir interaction energy. Within 
the worldline picture, it is already intuitively clear that for generic surfaces at a (suitably 
defined 2 ) distance a the Casimir interaction energy for a Dirichlet scalar is negative and a 
monotonously increasing function of a; therefore, the resulting force is always attractive in 
agreement with a recent theorem [3*Sj . 

2.2 Worldline numerics 

For the numerical evaluation of the expectation value Eq. (pj), two discretizations are 
required: first, the path integral is approximated by a finite sum over an ensemble of 
random paths x^(r), i = 1, . . . , til, each of them forming a closed loop in space(-time). 
Second, the propertime which parameterizes each path is discretized: 

Mr), re[0,T] Xik := xt(h ■ T/N), k = l,...,N; (12) 

2 A useful definition may be given by the following construction: let a > be the maximally possible 
distance between two auxiliary parallel plates that can be placed in between the surfaces constituting the 
Casimir configuration without mutual intersection. This excludes pathological cases such as surfaces which 
are folded into each other. In this construction, it should also be understood that a change of a should 
not be accompanied by a rotation of one of the surfaces. 
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i.e., the paths themselves are represented by N points per loop (ppl). Thus, the ensemble 
is described by a two dimensional array of space vectors (x^), with the indices i and k 
specifying the loop and the point on the loop, respectively. 

We generate the random paths using the v-loop algorithm j2Ej- This algorithm incor- 
porates the Gaussian term e _ 5 Jo dry ? as probability distribution, so that the path integral 
in Eq. (jSJ) becomes an arithmetic mean: 

/ Px(...)e-^o T ^ _ iVY..). (13) 

It is sufficient to generate only one so called unit-loop ensemble (yek), i.e., an ensemble of 
loops with center of mass xcm = and T = 1 . An ensemble with other values for xcm and 
T is then simply obtained by computing 

x^fc = xcm + Vf yek (14) 

for all I and k. At the same time, this technique provides for an analytic knowledge of the 
integrand's T dependence, which can be utilized for the T integration. 
With this discretization, the Casimir interaction energy Eq. (jlljl reads 

1 1 r°° dT f 1 x nL ^ 

^Casimir = --^^y 7^3 J <Px C M— ^ ©£ fccM + Vfy e ] . (15) 

The discretization error is controlled by the two parameters riL and N. The number of 
loops per ensemble «l is related to a statistical error of the arithmetic mean in Eq. (jlHj) . 
which can be determined by jack-knife analysis. The number of points per loop N is chosen 
sufficiently large to achieve the desired precision by studying the convergence of the result. 
In this work, we have used ensembles with up to = 2.5 • 10 5 and N = 4 ■ 10 6 . 

The major advantages of worldline numerics are its scalability and its independence 
of the background (geometry). The computational effort scales only linearly with the 
parameters nt,, N, D, etc, and the algorithm can be formulated for any given background 
geometry. A disadvantage is that the statistical error decreases only with l/y 7 ^, as for 
any Monte-Carlo method. This implies that high-precision computations may require high 
statistics, in contrast to estimates with, say, a few-percent error which require very little 
computational effort. 

For high-precision Casimir applications with an intended error of < 0.1%, the CPU 
time needed for the evaluation of Eq. (fT5|) can be reduced substantially by specializing the 
algorithm to the given Casimir geometry. Although this corresponds to a loss of generality, 
we believe that the strategies which we describe in the following, e.g., for the sphere-plate 
configuration, are examples for a general set of algorithmic tools which will be useful also 
for other Casimir configurations. 
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3 Worldline algorithms for Casimir configurations 



The general structure of a worldline algorithm for computing Casimir interaction energies 
is summarized by Eq. The only part of the algorithm that depends on the geometry 

consists of a diagnostic routine which checks whether a given loop (for given Xcm and T 
intersects with (more than one of) the surfaces S 1; S 2 , .... The result of this diagnostic 
routine immediately translates into the form of either the T or xcm integrand, depending 
on the actual order of integration. If the T integral is done first, the resulting loop-averaged 
xcm integrand can be viewed as the interaction energy density, the calculation of which is 
already an instructive intermediate step. 3 In principle, the T and xcm integration as well 
as the average over all worldlines can be done in arbitrary order, depending on numerical 
convenience. 

There is, however, an important technical difference between taking the worldline av- 
erage before or after the integrations. The apparent advantage of doing the worldline 
average first is that the resulting T and xcm integrands are smooth, despite the fact that 
the worldlines are fractal; this was exploited in many worldline numerical applications so 
far. In the present work, we nevertheless do the loop average at a later step. As a con- 
sequence, the resulting integrands can become complicated in the sense that the support 
of the integrand is a piecewise disconnected set. However, once the support is determined 
by special algorithms, at least one integral can be done analytically, since the integrand 
~ 0s = 1 and thus is extremely simple on the support. This leads to significant numerical 
acceleration, constituting the basic new ingredient of our improved algorithms. 

3.1 Sphere above plate 

The geometry of the sphere-plate configuration is illustrated in Fig. It is rotationally 
symmetric with respect to the z axis, therefore the three dimensional xcm integration in 
Eq. ((SJ) trivially reduces to a two dimensional integration. We choose the following order 
of remaining integrations/summations: first, we do the T integral for each worldline; then, 
we take the average over all loops, and finally integrate over the resulting energy density. 
For the first step, the numerically most challenging task is to determine the support Si of 
@s[ x cm + VTyp] on the T axis to perform the T integration in Eq. (fTo^l . In the given 
geometry, 0e[xcm + VTye] equals 1 if there exists a pair /c, /, such that x CM + VTyik lies 
inside the sphere and xcm+V^V^ lies below the plate; otherwise it is zero. 4 To investigate 
the support Si, it turns out to be useful to distinguish between xcm lying inside the sphere, 
x CM < R 2 , and xcm lying outside, x CM > R 2 , as the former case is much simpler than the 

3 The xcm integrand actually corresponds to a static effective-action density in the present case. The 
relation to the true interaction energy density which corresponds to the 00 component of the interaction 
energy- momentum tensor is give by a total derivative |32l I39j . 

4 Strictly speaking, this criterion misses the rare case that the link between two neighboring points which 
are both outside the sphere intersects the sphere. We neglect these contributions, since the verification of 
this pattern is much more time-consuming than simply increasing the amount of points per loop to reduce 
the corresponding systematic error. 
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Figure 1: Geometry of the sphere-plate configuration 



latter. 



3.1.1 Inside 

Inside the sphere, the support is a single T interval. The lower bound T™ in is given by 
the T value at which the loop xcm + VTye touches the plate, 

rr= ( R + a + x ™A \ (16) 

V mm fc y lkz ) 

where mm k yn kz is the minimal z coordinate of the unit loop yg. The upper bound T™ ax is 
the largest T value for which the loop intersects the sphere, 

Tr x = max{T : 3 fc (x CM + VTy ik ) 2 = i? 2 } (17) 



max 

k 



xcm • yek | / / x CM ■ y ek V x CM - R 2 



Performing the T integration, we obtain the Casimir interaction energy density ^casimir 
inside the sphere, £ C asimir = / d 3 x C M ^Casimir, 

1 1 ™ L / 1 1 \ 

^Casimir(xcM) = ^^2~ ^2 ^Jmaxj2 ~~ (T^ min ) 2 J ' ~ ^l^)' ( 19 ) 

where the # function takes care of the (non-generic) case that the loop never intersects 
both surfaces. This quantity is plotted in the contour plot Fig. El in the region inside the 
white circle. The contribution to the total Casimir interaction energy is small compared 
to the energy outside the sphere. The density of the latter is shown in the same figure 
outside the circle, obtained by the procedure described next. 
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Figure 2: Contour plot of the negative Casimir interaction energy density ecasimir for a 
sphere of radius R above an infinite plate; the sphere-plate separation a has been chosen 
as a = R here. 

3.1.2 Outside 

Outside the sphere the support Se is not merely one single interval as in the previous case, 
but a whole set of successive intervals. As illustrated in Fig. El for a unit loop ye at a center 
of mass x CM outside the sphere, the ray x CM + VTyik does not pierce the sphere for most 
indices k. The corresponding points on the loop are not relevant for the Casimir energy 
and the first step in our algorithm is to sort them out. Two conditions are evaluated for 
this purpose: a point yek is only relevant for further computations if 

1. the vector yek points towards the sphere, implying 

x C m • yek < 0, (20) 

2. the distance h between the ray x CM + VTyek and the center of the sphere is smaller 
than the radius R, 

h 2 = x CM - ( x CM • < R 2 . (21) 

V \yek\J 

If these conditions are fulfilled, the T values at which the ray intersects the sphere are 
determined by 

(x CM + Vf y ek ) 2 = R 2 , (22) 
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Figure 3: For any propertime T, the ray xcm + VTytk does not intersect the sphere. The 
corresponding point y# is thus not relevant for the interaction energy density at the given 
xcm and consequently sorted out in a first step of the algorithm. 



which has the solutions 




^H-=^±lP^ --"^T^ ■ (23) 



For T e [T^, T^j~] the point x CM + \fTy ik lies inside the sphere and consequently we know 
that the loop intersects the sphere. For a given loop, the total set of T values for which 
this is the case is the union of the intervals \T[ k) T^ k ] of all points in the unit loop, 
Ufct^ffci^fc]- Taking into account the minimal T value for which the loop intersects the 
plate, T™ m , the contribution of the unit loop to the propertime integrand has the support 

Sz=[Tf n ,oo)n\J[T ik ,T i t}. (24) 

k 

The set union can be determined efficiently by use of a sorting algorithm like quicksort, 
for example. Once Si is determined, the T integration can be performed analytically. The 
worldline estimate for the Casimir interaction energy density outside the sphere therewith 
is 

£Casimir(x C M) = J s ^3' ( 25 ) 

which is plotted in Fig. |21 outside the white circle. 



3.1.3 Optimization 

The algorithm so far works well if the distance between sphere and plate a is of the 
same order of magnitude as the sphere's radius R, a w R. To improve the accuracy by 
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increasing the number of loops and the number of points per loop N, the algorithm 
can be parallelized as embarrassingly parallel computation by dividing the loop ensemble 
into independently processed sub-ensembles. However, if the two scales a and R differ 
significantly in size, additional improvements of our algorithm are advisable. 

Large distances (a ^> R): if the distance a is large compared to the radius R, the algorithm 
described so far becomes inefficient due to the following reason: only loops with a minimal 
extent of the same order of magnitude as the distance between sphere and plate do con- 
tribute to the Casimir energy. For a loop x CM + VTye this means, that the unit loop yt 
has to be scaled by a large factor y/T. This implies that the distance between subsequent 
points on the loop increases, too. However, to ensure that the scaled loops still resolve the 
sphere, this distance should be significantly smaller than the sphere's radius. Thus, the 
number of points per loop N has to be increased with increasing a/R. 

A rough measure for the extent of a loop is the variance of the coordinates of its points. 
The ensemble average of the variance for large N is ((VTy ek ) 2 ) = T/6. As a consequence 
we expect the T integral to be dominated by T « 6a 2 , also because the contribution for 
large T is damped by the 1/T 3 factor. The root-mean-square of the distance between two 
subsequent points on a loop for large N is a = \J2T/N . Using the dominating T value 
we obtain a « 2\j3/Na. We demand this value to be much smaller than the radius of the 
sphere, which implies N ^> 12a 2 /R 2 . For a distance a = 10R, already much more than 
1000 ppl have to be used, for a = 100R much more than 100 000 ppl. 

A slight modification enables our algorithm to cope with this high resolution and the 
corresponding amount of data much more efficiently. So far, for all center of masses xcm 
with a common z coordinate Xcm z , the first interval on the right-hand side of Eq. (|24j) . 
[T™ m , oo), is the same, which can be utilized to speed up the calculation. In contrast, the 
union in the same equation is different for all centers of masses. However, it is this part 
of the equation which consumes most of the CPU time. Modifying the transformation Eq. 
dHj) reverses the circumstances: let us define the rotation i?(x CM ,e 2 ) by x CM /|x CM | = 
i?(x CM , e z )e z . By using 

x £fc = x CM + ^(xcm, e z )ym (26) 

(see Fig. HJ), the union in Eq. (}2"4"j) is the same for all center-of-mass values xqm with 
a given absolute value |xcm| and can be computed once and for all. In turn, T™ m is no 
longer degenerate with respect to some xcm coordinate. The important advantage is that 
this dependence can be computed much faster. For each loop, we generate an array of 
its minimal z coordinate as function of the angle between xcm and e z . The bound T™ m 
then results from Eq. (|16|). where the minimum is read from the array. Note that the 
transformation f|26|) is a legitimate symmetry operation for ensemble- averaged quantities, 
owing to the rotational invariance of the exponential weight factor in Eq. (|13p. 

There is a price to be paid for the desired feature of having a common set union in 
Eq. (J24)) for different centers of masses £cm : without the transformation, all points of 
a given unit loop are equally involved in scanning the curvature of the sphere. With the 
transformation, always the same points of a unit loop are close to the sphere, independently 
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Figure 4: Sketch of a unit loop at different centers of masses xcmi and xcm2, as used for 
large a/R. The unit loop is rotated corresponding to the orientation of the center of mass. 



of xcm- This corresponds to a loss of statistics, implying a slight increase of the statistical 
errors. However, this is by far compensated for by the gain in computation speed, which 
enables us to significantly reduce again the statistical error by brute force. 

Small distances (a ^ R): the main contribution of the Casimir interaction energy density 
is localized between sphere and plate. If the distance a is much smaller than the sphere's 
radius R, the lower bound of the support Se in that region is at very small T values 
compared to the upper bound of the support's first interval. Since the T integrand falls 
off rapidly with 1/T 3 , the T integral is dominated by this lower bound. Therefore, a very 
good estimate is given by replacing Se simply by the interval [T™ m , oo), resulting in 



"Casimir (XCM) ^ "g^^ ( T min)2 " ^ 

£=1 ^ 



In particular outside the sphere, the numerical evaluation of this expression is much faster 
than the evaluation of the full expression Eq. (|25|) . 

For a rough estimate of the validity range of this approximation, we use the ensemble's 
standard deviation of the point position, ((VTyek) 2 ) = T/6, to estimate the extent of 
a loop. Right between sphere and plate, where the energy density is largest, the lower 
bound of the propertime integral then is approximately T™ m ps 3a 2 /2. If the extent of the 
loop increases beyond a/2 + 2R, we expect the loop to intersect the sphere no longer for 
T\ <; 6 (a/2 + 2R) 2 . By setting the value of T\ to infinity instead, as done in Eq. (|27j). we 
introduce a systematic error Aecasimir/^Casimir ~ (l +4R/a) 4 . For distances a smaller 
than 0.8R this error is smaller than one per mille. In this work, we have used Eq. (127)) to 
compute high-precision values for a < 0.02R. 
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Figure 5: Casimir interaction energy of a sphere with radius R and an infinite plate vs. 
the curvature parameter a/R. The energy is normalized to the zeroth-order PFA formula 
(J2J), Epp A . For larger curvature parameter, the PFA estimate (dot-dashed line) differs 
qualitatively from the worldline result (crosses with error bars). Here, we observe good 
agreement of our result with the exact solution of j2Hj which is available for a/R > 0.1 
(dashed line). 



3.1.4 Results 

Figure 03 presents a global view on the Casimir interaction energy for a wide range of the 
curvature parameter a/R; the energy is normalized to the zeroth order of the PFA formula 
(J2J), Epp A . For small a/R ("large spheres"), our worldline result (crosses with error bars) 
and the full sphere- and plate-based PFA estimates (dashed-dotted lines) show reasonable 
agreement, settling at the zeroth-order PFA E p °p A . The full PFA departs on the percent 
level from E PPA for a/R > 0.01, exhibiting a relative energy decrease. By contrast, our 
worldline result first stays close to E PPA and then increases towards larger energy values 
relative to E p °p A . This observation confirms earlier worldline studies j2Ej and agrees with 
the optical approximation [27] in this curvature regime. 

For larger curvature a/R > 0.1 ("smaller spheres"), we observe a strong increase relative 
toE P % A [33]. Here, our data satisfactorily agrees with the exact solution found recently for 
this regime [2E] (dashed line). The latter work also provides for an exact asymptotic limit 
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Figure 6: Magnified view of Fig. for small a/R. The 0.1% validity range of the PFA 
is characterized by curvature parameters, where the error band of our worldline results 
and the PFA band (blue-shaded/in between the dot-dashed lines) overlap, see Eq. (|28|). 
The dashed line depicts a constraint polynomial fit of the worldline result, p(a/R) = 
1 + 0.35 (a/R) — 1.92(a/R) 2 , and its standard deviation, see Eq. (|3U|). The inlay displays 
the same curves with a linear a/R axis. 



for a/R — > oo, resulting in 180/7T 4 for our normalization. Our worldline data confirms this 
limit in Fig. 

Two important lessons can be learned from this plot: first, the PFA already fails to 
predict the correct sign of the curvature effects beyond zeroth order, see also |3D] . Second, 
the relation between the Casimir effect for Dirichlet scalars and that for the EM field is 
strongly geometry dependent. For the parallel-plate case, Casimir forces only differ by the 
number of degrees of freedom, cf. the coefficient cpp in Eq. (|TJ. For large curvature, the 
Casimir energy for the Dirichlet scalar scales with a -2 , whereas that for the EM field obeys 
the Casimir- Polder law ~ a -4 |U E2] • Already this difference demonstrates that simple 
approximation methods such as the PFA are highly problematic, since no reference to the 
nature of the fluctuating field other than the coefficient cpp is made. 

For a quantitative determination of the PFA validity limits, Fig. El displays the zeroth- 
order normalized energy for small curvature parameter a/R. Here, our result has an 
accuracy of 0.1% (jack-knife analysis). The error is dominated by the Monte Carlo sampling 
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and the ordinary-integration accuracy; the error from the worldline discretization is found 
negligible in this regime, implying a sufficient proximity to the continuum limit. 

In addition to our numerical error band, we consider the region between the sphere- 
and the plate-based PEA as the PFA error band. We identify the 0.1% accuracy limit of 
the PFA with the curvature parameter a/R\ ,i% where the two bands do no longer overlap. 
We obtain 

a PFA 

- < 0.00073 (28) 
R o.i% 

as the corresponding validity range for the curvature parameter. For instance, for a typical 
sphere with R = 200/im and an experimental accuracy goal of 0.1%, the PFA should not 
be used for a > 150nm. We conclude that the PFA should be dropped from the analysis 
of future experiments. 

For the 1% accuracy limit of the PFA, we increase the band of our worldline estimate 
by this size and again determine the curvature parameter for which there is no intersection 
with the PFA band anymore. We obtain 

a 
R 

For a sphere with R = 200/um and an experimental accuracy goal of 1%, the PFA holds 
for a < 1.5/xm. This result confirms the use of the PFA for the data analysis of the 
corresponding experiments performed so far. 

In order to study the asymptotic expansion of the normalized energy, we fit our worldline 
numerical data to a second-order polynomial for a/R < 0.1 and include the exactly known 
result for a/R — > 0. We obtain 



PFA 

< 0.00755. (29) 

i% 



«t = -rr^f ( 1 + 0.35^ - ± 0.19^1 - 137.2^ + 5125^ 

(30) 

valid for a/R < 0.1; here, cpp = 1 for the real and cpp = 2 for a complex Dirichlet scalar. 
The fit result is plotted in Fig. |H1 (dashed lines), which illustrates that E ~ £wn data fit 
is a satisfactory approximation to the Casimir energy for a/R < 0.1, replacing the PFA 
The inlay of Figure El displays the same curves with a linear a/R axis, illustrating 
that the lowest-order curvature effect is linear in a/R. A more direct result for the linear 
curvature coefficient can be obtained by a constraint linear fit; in this simpler case, the fit 
polynomial yields p^t{a/R) = 1 + (0.33 ± 0.06) ^ instead of the expression in parentheses in 
Eq. (jHOJ). Given the results of the PFA (JTJ), the semiclassical approximation [21], p sc (a/i?) ~ 
1 — 0.17-|, cf. [2E], an d the optical approximation p op t(a/R) — 1 + 0.05-|, the latter 
appears to estimate curvature effects more appropriately; but all these approximations are 
not quantitatively reliable for beyond-zeroth-order curvature effects. 

3.2 Cylinder above plate 

The cylinder-plate configuration is a promising tool for high-precision experiments |43j . 
since the force signal increases linearly with the cylinder length. The numerics is very 
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Figure 7: Casimir interaction energy (normalized to -Epfa) °f an infinitely long cylinder 
with radius R at a distance a above an infinite plate vs. the curvature parameter a/R. We 
observe good agreement of our result with the exact solution of [22] which is available for 
a/R > 0.1 (dashed line). 



similar to the sphere-plate configuration, even less computing power is required, because 
only two dimensional loops have to be processed due to the translational symmetry. Figure 
13 shows the corresponding Casimir interaction energy versus the curvature parameter. The 
energy axis is again normalized to the zeroth-order PFA result, 

(31) 

A magnified view of the small curvature region in Fig. [7| is shown in Fig. |HJ As for the 
sphere-plate configuration, we fit our data to a second-order polynomial in this range, 
including the exactly known result for a/R — > 0, yielding 



E??* n tl?^ e = --^= ^ [l + 0.2l4 - 0.66^ ± 0.0974a /l - 68.60^ + 1282° 



1920\/2 a 5 / 2 V" ' *""*R R 2 ^ """" ' R V * """"" R ' ^ R 2 

(32) 

for a/R < 0.1. The inlay of Fig. El shows the same data with a linear a/R axis. As for 
the sphere-plate geometry, this plot demonstrates that the lowest-order curvature effect is 
linear in a/R. A simpler linear fit to our data results in p&t(x) ~ 1 + (0.195 ± 0.028)-^. 
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Figure 8: Magnified view of Fig. [7| for small values of a/R. The dashed line with error 
band depicts a constraint polynomial fit to the numerical data, p(a/R) = 1 + 0.21 (a/R) — 
0.66(a/R) 2 , and its standard deviation. The inlay displays the same curves with a linear 
a/R axis. 



This is in remarkable agreement with the recently found analytical result p(a/R) = 1 + 
0.194|r + 0(a 2 /R 2 ) [20]; which represents a strong confirmation for both methods. 

The qualitative conclusions for the validity of the PFA are similar to that for the sphere 
above a plate: beyond leading order, the PFA even predicts the wrong sign of the curvature 
effects. Quantitatively, the PFA validity limits are a factor ~ 3 larger than Eqs. (|28 j) . (|29| ) . 
owing to the absence of curvature along the cylinder axis. 

The most important difference to the sphere-plate case arises for large a/R. Here, the 
data is compatible with a log-like increase relative to E p °p A , implying a surprisingly weak 
decrease of the Casimir force for large curvature a/R — > oo. Our result agrees nicely with 
the recent exact result [2H] which is available for a/R > 0.1. The data thus confirms 
the observation of (2H] that the resulting Casimir force has the weakest possible decay, 
F ~ l/[a 3 \n(a/ R)], for asymptotically large curvature parameter a/R — > oo. 
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Figure 9: Geometry of Casimir's parallel-plate configuration. A loop contributes to the 
Casimir interaction energy if its extent L(_ along the z direction is larger than the plate 
distance a. 

3.3 Parallel plates revisited 

As discussed at the beginning of this section, the order of the T and xcm integration and 
the ensemble averaging can be chosen arbitrarily. As an example for an "unusual" order, 
let us reconsider Casimir's classic parallel-plate configuration in D = d + 1 dimensional 
spacetime, doing the xcm integral first and keeping the ensemble average till the very end. 
This will reveal an unexpected mapping between the .D-dimensional Casimir effect and 
standard polymer physics. 

In d space dimensions, the surface or area volume A of the Casimir plates is taken 
as d — 1 dimensional. The two (hyper-)plates are separated by a distance a along the z 
direction which is normal to the plates, see Fig. El For this configuration, the Casimir 
interaction energy for the massless Dirichlet scalar boils down to 

^ ftL POO ^rp poo 

^Casimir = ~A 2 ^ D/2 J Q J^UJl J ^ ^ CM 6 ^ CM + ' ( 33 ) 

where y Z £ denotes the z coordinate of the £th unit loop. Let us denote the extent of the 
£th unit loop in the z direction by Li, 

L e := max (\y zU - y ztk \) , (34) 

i, k 

see Fig. El A scaled unit loop intersects both plates if VTLg > a. For a given unit loop with 
extent L e and for a given propertime value T, the support of the zqm integral corresponds 
to an interval I Z (T, Lg) = s/TLt — a. Independently of the precise location of this interval 
on the zcm axis, the ^cm integral yields, 

/oo 
dz CM Q[z C m + VTy z£ ] = (VTL £ - a) 6(VTL e - a). (35) 
-oo 
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Now, also the T integral can be done analytically, resulting in 



A 



Casimir 



D 



a " 1 D(D — l)(47r) D / 2 
A 1 



<L?>. (36) 



a D-i D(D - l){4n) D /< 



We observe that the Casimir interaction energy of the parallel-plate configuration in D = 
d + 1 spacetime dimensions is proportional to the -Dth moment of the ensemble-averaged 
extent of a unit loop. This ensemble average could easily be performed, which would lead 
us back to the results of [2S|. 

Here, we will be satisfied by hilighting Eq. (}36|) from a different perspective, namely, 
in the language of polymer physics. The GauBian velocity distribution of our worldlines 
is identical to the Hamiltonian of a polymer, i.e., the continuum limit of a random chain, 
without self-avoidance or excluded-volume effect jH] in the limit of zero end-to-end dis- 
tance. In this language, Li corresponds to the maximum spatial extent of the closed 
polymer measured in units of yl v c/ (2D p ). Here, l p = Nc denotes the total length of the 
polymer, c is the chain length, and D p is the number of dimensions in which the polymer 
can move; the latter is completely arbitrary and independent of the dimensionality of the 
Casimir system. Since Le is a highly nonlocal object, its ensemble average is actually not 
so easily computable by standard methods. Our result Eq. ()36|) now maps the problem 
of computing any Dth moment of Li on the D- dimensional Casimir problem for parallel 
plates. Using the standard result for the latter [35| . 

Casimir = ~ jjy^ T(D/2) ((D), (37) 

we obtain by comparison with Eq. (J36)) . 

(L?[y]) = D(D-l)T(D/2)((D), (38) 

a result that we have so far not been able to find in the literature of polymer physics. 
Even the limit D — > 1 can be taken, corresponding to the Casimir effect in zero space 
dimensions: this results in (L([y\) = y/n for the average extent of a closed polymer. 



4 Conclusions 

We have presented improved worldline numerical algorithms that can efficiently deal with 
Casimir configurations involving curved surfaces. We have used these algorithms to com- 
pute Casimir interaction energies for the sphere-plate and cylinder-plate configuration in- 
duced by a scalar field with Dirichlet boundary conditions. These computations are done 
from first principles for a wide range of curvature parameters a/R. In general, we ob- 
serve that curvature effects and geometry dependencies are intriguingly rich, implying that 
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naive estimates can easily be misguiding. In particular, predictions based on the PEA are 
only reliable in the asymptotic no-curvature limit with quantitative validity bounds given 
above. We have constructed polynomial fits of our results which can be used in the small- 
curvature regime, a/R < 0.1, as a well-founded substitute for the PEA formulas. Given 
the size of the true curvature corrections for the Dirichlet scalar, we expect that genuine 
Casimir curvature effects are in reach of currently planned experiments. In this spirit, 
the so-called lateral Casimir force for corrugated surfaces has recently been proposed as a 
suitable candidate for identifying non-trivial geometry dependences beyond the PEA |46j . 

Beyond the Dirichlet scalar investigated here, it is well possible, e.g., for the EM field, 
that some cancellation of curvature effects occurs between modes obeying different bound- 
ary conditions. In fact, such a partial cancellation between TE and TM modes of the 
separable cylinder-plate geometry can be observed in the recent exact result for the EM 
field for medium curvature for small curvature, curvature effects can even reverse 
sign j3D]. More quantitatively, the TM mode in the cylinder-plate case obeys Dirichlet 
boundary condition and thus contributes, e.g., to the linear curvature correction with a 
coefficient ~ 0.194, as discussed below Eq. (|52J): the TE mode obeys Neumann bound- 
ary conditions, giving a negative contribution which in total turns this linear coefficient 
for the EM field into ~ —0.48 [3U]. The latter result, in fact, lies in the broad range of 
[—0.92, —0.25] spanned by the PFA; since the PFA does not make any reference to the 
nature of the fluctuating field, this rough coincidence is, of course, purely accidental. This 
strong dependence of Casimir curvature effects on the nature of the fluctuating fields alone 
demonstrates already that approximations ignoring this difference such as the PFA cannot 
be trusted. We emphasize again that Casimir calculations for the EM field in non-separable 
geometries, such as the important sphere-plate case, remain a prominent open problem. 

From a technical point of view, we would like to stress that our results demonstrate the 
capability of worldline numerics for performing high-precision computations with compar- 
atively little computing power. The simple scalability of the algorithms and the flexibility 
for adapting them to arbitrary geometries makes worldline numerics a unique tool for 
computing quantum energies. 

Our algorithmic strategies also revealed an unexpected mapping between the D-dimen- 
sional parallel-plate Casimir effect and aspects of a random-chain polymer ensemble. The 
origin of this mapping, of course, lies in the fact that both quantum fluctuations in Casimir 
systems as well as a polymer ensemble can be described by Feynman path integrals. In 
the present case, the mapping can be utilized to transform a comparatively difficult poly- 
mer problem into a field-theoretic Casimir problem which can be solved by a variety of 
techniques. We believe that this mapping is just a special case of a more general class of 
mappings with potentially fruitful applications in both directions. 
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